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Abstract 

Motivated by the goal to give the simplest possible microscopic foundation for a broad 
class of topological phases, we study quantum mechanical lattice models where the 
topology of the lattice is one of the dynamical variables. However, a fluctuating ge- 
ometry can remove the separation between the system size and the range of local in- 
teractions, which is important for topological protection and ultimately the stability 
of a topological phase. In particular, it can open the door to a pathology, which has 
been studied in the context of quantum gravity and goes by the name of 'baby uni- 
verse'. Here we discuss three distinct approaches to suppressing these pathological 
fluctuations. We complement this discussion by applying Cheeger's theory relating the 
geometry of manifolds to their vibrational modes to study the spectra of Hamiltonians. 
In particular, we present a detailed study of the statistical properties of loop gas and 
string net models on fluctuating lattices, both analytically and numerically. 
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1. Introduction 

Topological states of matter have become a major topic in condensed matter theory 
and experiment. Quite simple microscopic Hamiltonians, e.g. those of fractional quan- 
tum Hall (FQH) systems, are known to harbor emergent topological phases but often 
determining the correct effective low-energy theory is diflicult: it becomes a delicate 
matter of energetics once one moves beyond idealized short ranged interactions. An- 
other way of idealizing interactions comes in the form of quantum lattice models. In 
recent years lattice models have become another, complementary source of topological 
phases. Examples include the well-studied toric code model yj] and the Kitaev hon- 
eycomb model [2], which realizes the Ising topological quantum field theory (TQFT) 
in a controlled perturbative regime. A more general class of lattice models was intro- 
duced by Levin and Wen ^ - "Turaev-Viro" 101 in the math literature - which are rich 
sources of exactly solved achiral ("doubled") theories. 

There is, however, a trade-off' in using these models. In the exactly solved LW 
lattice models we know everything, but the Hamiltonians look artificial: finely tuned, 
I2-body interactions without the comfort of familiar "kinetic" and "potential" terms. 
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In FQH the Hamiltonian is simple and natural, but the connection to the (or several 
competing) topological phases is not obvious. The motivation of this paper is the search 
for a middle ground where a simple local Hamiltonian can be clearly identified with a 
unique topological phase. Our chief innovation is to treat the lattice itself as a dynamic 
variable so that there is no fixed lattice underlying the Hilbert space, hence the words 
"off lattice" in the title. Our search for this middle ground is presented as a travel 
log, with some surprises, disappointments, and discoveries. Along the way we came 
to better appreciate what exactly a fixed lattice is good for and what adaptations its 
absence requires. Briefly, a lattice model supplies two length scales, the lattice scale a 
and the length of period L, where topological protection comes from an error scaling: 
e ~ g-const-L/ci Protecting quantum information without this ratio of scales, e.g. when 
going off lattice, is a key challenge. 

Levin- Wen models place degrees of freedom ("labels") on the edges of a trivalent 
graph dual to a fixed triangulation A of a surface Z, where the labels come from a 
fusion category. The labeled graph is called a "string net". There is a mathematical fact 
encouraging us to leave the concept of a fixed lattice behind. Starting with a consistent 
F symbol (i.e. one obeying the pentagon relations) for a fusion category 



(1.1) 



and a closed surface 2, the vector space of admissibly labeled string nets on 2 modulo 
isotopy and F is canonically isomorphic to the LW ground state Hilbert space V(Z) built 
from a fixed dual triangulation. Clearly ( 11.1) is a simpler starting point than the LW 
Hamiltonian, but there is a serious problem: no single lattice (we use the word "lattice" 
here to mean the dual structure to a triangulation of E) will contain simultaneously 
bonds forming the "stick figures" on both the left and right hand sides of (11.11 ): the 
F-move necessarily modifies the lattice. 

In fact the LW Hamiltonian arises from a composition of such F-moves which — 
like a perturbation through excited states — begins and finally ends on the same lattice. 
For the honeycomb lattice, the LW Hamiltonian is a 12-body term, arising from 6 F- 
moves. Our idea here is to construct an environment where the F-move itself may be 
written directly into the Hamiltonian forcing the underlying lattice to fluctuate. 

This paper starts by reviewing the problems we encounter when going off-lattice 
in Sec. |2]and then present various ways to achieve off-lattice topological models in 
Sec. [3]before concluding with a summary of the lessons learned from this off-lattice 
expedition in Sec. E] The appendices contain a number of results that are interest- 





ing beyond the main subject of this paper Appendix A proves a discrete version of 



Cheeger's theorem appropriate graphs with weights on edges and vertices extending 
work of Chung [50. We use it to derive bounds on the gaps of graph Laplacians and 



other local Hamiltonians. Appendix B derives an upper bound for the scaling of the 



gap of the graph Laplacian on outer planar triangulations, following an idea developed 



with Oded Schramm. Appendix C presents exact results for an off-lattice version of 



the toric code model, and finally Appendix D presents numerical results for the gap of 



the Graph Laplacian and Cheeger's constant for trees, triangulations of the sphere and 
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outer planar triangulations. 



2. Off Lattice Hazards 

Let us begin by cataloging three hazards which await us off-lattice. 

2.1. Baby Universe 

The first hazard is well known in Euclidean quantum gravity. It is called "minbus" 
or "baby universe" [6]. It refers to the fact that if a triangulation (or any other poly- 
hedral decomposition) is chosen uniformly at randorcQ for a sphere (or other closed 
surface of fixed genus g) among all triangulations of a fixed number n of triangles, it 
is likely that there will be a short dual loop containing numerous vertices on both sides 
(or in the case of surfaces with g > there will also be short nonseparating loops). We 
measure the length of a (dual) loop simply as the number of edges it crosses. This can 
be formalized by saying that Cheeger's |7] isoperimetric constant — > when n — > oo, 
almost always 

min ^^"^^'^^^ _ , (2.1) 

separating dual loops y mm[area{S), area{S)] 

where S and S are the components of the surface E minus y, and area(S) (area(S )) is 
the number of vertices (or sites) in S (S). 

It is known that typically a triangulated sphere has k x (9(^^). There is in fact an 
asymptotic formula [8] for the number #g(n) of isomorphically distinct divisions of a 
closed genus g surface into squares with the condition that the dual graph has no odd 
cycles: 

#g(n) = 12"ni'-«-^^ = 12"n"5^, (2.2) 

where x is the Euler number It is only a simplifying technicality to treat these quad- 
rangulations rather than triangulations — similar asymptotics should apply but with 12 
replaced by some other less convenient base; the n^^^ is universal. 

An easy application of Eq. ( 12.2b (we thank Gilles Schaeffer for bringing this to our 
attention) is that it is possible to estimate the fraction /4(n) of genus g surfaces that are 
divided by a separating curve of length four into surfaces of genus and g2, gi,g2 > 0. 
The method is to remove a random square from random surfaces of and g2, glue the 
results together and count how many ways this is possible as a fraction of surfaces of 
genus g - gi + g2- The result is: 

/4(n) « n-'^\ (2.3) 

To summarize: constant size bottle necks are algebraically likely and logarithmic 
bottlenecks are virtually assured. This appears to be very unfavorable for the protection 
of topological information. When working with lattice models, we are used to error 
rates appropriate to tunneling problems like e^'^°"'^^'^, where L is the linear dimension 
of the lattice as a multiple of the lattice constant. If bottlenecks reduce L to constant or 



or by any other local formula. 
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even logarithmic size, the protection disappears. This is the first issue to come to terms 
with when considering uniformly random triangluations (URT). 

For quite a different reason, these bottlenecks and baby universes have been unwel- 
come also in quantum gravity. There one seeks a Hamiltonian function on triangula- 
tion which concentrates near Euclidean flat space but still allows a liquid of possible 
universes. Forty years of effort have failed to find such a phase even in dimension 
2 (at least in the homogeneous setting), for a review see It is possible to form 
a "branching polymer universe" of Hausdorff dimension 2, to perturb about a single 
rigid Euclidean crystal, or a "collapsed phase" of infinite Hausdorff dimension, but a 
nearly flat, yet liquid, phase, or even a critical point, has been elusive. However, in the 
last dozen years progress has been made by breaking the symmetry between space and 
time and allowing only triangulations appropriately foliated by space-like leaves |@]. 
The approach is called "causal dynamical triangulation" (CDT) and has been shown 
numerically to provide "birth control" [|9J — there are parameter regimes, called the 
"C phase", with no baby universes in which the space-like leaves are on average nearly 
Euclidean of the desired dimension. In Section 2.3 we describe an approach to build- 
ing a(l-i-l) + l dimensional model for a 2+1 -dimensional anyonic system in which 
the 2 spatial dimensions are broken into a 1 + 1 pair to exploit the favorable statistical 
geometry of 1 + 1 CDTs. 

2.2. Gapless Modes 

The second issue with URT is the mixing time. We analytically estimated the 
Cheeger isoperimetric constant h in a toy model of surface triangulations called "outer 
planar" triangulations. In this context we show analytically h ^ yTT-, whereas in our nu- 



merical study (presented in Appendix D the first eigenvalue A of the graph Laplacian 
(the "graph" has vertices outer plane triangulations and edges plaquette flips between 
these) goes like A ^ rc'^ (this translates to the first eigenvalue A of the graph incidence 
matrix scaling like A a since A k An^^ in our model^. Cheeger like inequalities 



show (see Appendix Ai: 



2/z ^ /i ^ y , (2.4) 

or 

)p h )p yT'^ . 

We believe the truth is near the high end h ^ n and that on the sphere A ^ n 

As a further probe of the spectrum, we studied the dynamics of string nets on the 
2-sphere (see Appendix D i and also observed a mixing time n^, i.e. A x n^^, and 
A X n String nets are dual to triangulations but slightly more flexible, e. g. a closed 
loop in a string net is permitted whereas the usual definition of triangulation does not 
permit a triangle to be glued to itself. 

All this confirms the findings of the quantum gravity community: the space of 
random triangulations, quadrangulations, string nets, etc. on a surface will mix alge- 
braically fast but not so fast (which would need to be 0{n)) so that the A first eigenvalue 



^ We need to remove an exponentially small number of states from the Hilbert space to obtain this scaling. 
For a detailed discussion see |Appendix C| and |Appendix D| 
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of the incidence matrix is gapped or, equivalently, that the Cheeger constant k (of the 
graph-of-triangulations) is bounded away from zero. 



2.3. Local Distinguishability 

A third problem was noticed when we studied muhi loops on a periodic honey- 
comb lattice. The chosen dynamics is that of the toric code |1], however, we made 
the convention that if two multi loops were isotopic, deformable one to the other, then 
they would be identified and represented by a single ket. In Appendix C we present 
data which shows that the trivial winding sector can easily be picked out from the other 
three by a "local" observation — we count the number of "leaves" — that is loops with 
no smaller loops within. In retrospect, this is no surprise. Being in the trivial sector 
allows the possibility of no essential loops whatever — this possibility permits more 
space on the lattice for leaves. 

Local distinguishability is, of course, the death knell of topological protection. A 
state which can be observed locally can be acted on by a local operator. This is a third 
disturbing finding if we grant that leaves are to be considered local structures. Since 
metrical notions have been temporarily banished, it is up to our intuition to reformulate 
the appropriate meaning of "local opeator." Leaf detection has, in this context, as good 
of a claim to being local as does any operator. 



2.4. Work-arounds for off Lattice Troubles 

Of the three problems, (2) is the least concerning. Even if there are low energy 
metrical fluctuations (one may dub them "gravity waves"), they appear decoupled from 
topological degrees of freedom which can be encoded on each lattice. With respect to 
problems (1) and (3), there is a somewhat solipsistic solution to the apparent loss of 
topological protection from (1) bottlenecks and (3) variations in leaf count. It is simply 
to deny that this is a problem. Once kets of the Hilbert space are isotopy classes (of 
triangulations or nets — perhaps together with a particle type labelling of the bonds) 
we have lost direct contact with any notion of a position coordinate x. Isotopy slides 
and stretches, so we no longer know what is long or short or even where we are. This 
viewpoint leads one to say there are no local operators at all and therefore topologi- 
cal protection — protection against local operations — is tautological. But such a view 
comes with a heavy price — without a position coordinate x correlation functions loose 
meaning and contact with the condensed matter notions disappears. Consequently we 
will not take this path but rather consider three distinct approaches all of which en- 
force flat Euclidean space as the background, but in quite different ways and with quite 
different results. 

In summary, we find that there are plausible and even intriguing ways to model 
topological phases off lattice. The next step should be to identify a case where the 
model variables can be mapped to electron degrees of freedom. 



3. Enforcing Flat Space: Crystalline, Liquid off Lattice Models, and CDT 

Topological phases of lattice models are known in the physics literature from Levin 
and Wen [3] and in the math literature from Turaev and Viro [4]. We explore what 
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happens in taking such a model "off lattice" by including the underlying cell structure 
or "lattice" among the dynamic variables. Thus, our Hilbert space 'H will be spanned 
by kets which are pairs |(A, 5)) where A is a triangulation of the surface S with n 
(fixed) triangles and 5 is a "labeling" of the dual edges of A. We can equally well 
focus on the dual string net and its variations, some of which may not be dual to actual 
triangulations. It is the dual edges which form our so-called lattice, e.g. the dual edges 
form a honeycomb if A is the standard triangulation of the plane by equilateral triangles. 
The labels are from a pivotal fusion category. Two interesting examples are with label 
set { 1 , r) of the Fibonacci theory Fib and fusion rule t®t — 1 © r, which when applied 
as labels on string nets yields the theory Dfib lllOll . or with label set {l,x) and fusion 
rule X® X — 1, which gives the toric code f^. Rather than speaking in generalities, 
we give our constructions in the former case. They are easily extended to the broader 
class. 

In moving off lattice we can be timid or bold but as we have argued we must find 
some way to tie our lattices to the Euclidian plane. Our Hamiltonian can charge energy 
for defects in a base lattice, say the honeycomb, or it can treat all triangulations equally. 
There is, of course, an adjustable parameter connecting large to zero energy penalty 
for defects. However, decades of experience with 2D quantum gravity (qg) models 
suggest that there is a single phase transition from a phase with amplitudes clustered 
near the original honeycomb (we call this crystalline) and a liquid of lattices whose 
geometry is almost surely "cactus shaped" - Cheeger constant ~ jn the quantum 
gravity community the lack of an intermediated phase was the cause of some despair, 
the cactus buds being called "baby universes." As we remarked in the introduction, this 
problem may have been solved in the quantum gravity context [9, 6] by introducing an 
appropriate causal structure. This is explored in Section 2.3, where 3D space is split 
into a radial "space" coordinate p and a periodic "pseudo-time" coordinate 9. But if one 
insists that 2D space be treated homogenously, then the baby universes must be faced. 
We will do this but first let us explore moving only timidly off the honeycomb. There 
we find a gapped model which is conceptually very simple (we think more simple 
than the LW model) but the price of our timidness is that the gap is absurdly small, 

perhaps ~ ei-^j , where e is the energy scale of the individual terms. In spite of the 
disappointingly small gap, we next explain this model as it is a nice, controlled context 
for stepping - ever so slightly - off lattice. A bolder step will be taken later. 

3.1. Crystalline Case 

Our first Hamiltonian has the form Hgg - H^^ + 5D. Hqg acts on the direct sum of 
fibers of a bundle of states over the moduli space of string nets (thought of as metrics) 
on S, say a torus. The terms of H^^ are of two types: 

1. Fusion constraints; these are projectors acting within fibers 

2. F-moves; these act between adjacent fibers 

The fibers are degrees of freedom on the edge set of any given string net. F-moves de- 
fine a connection linking these fibers together which, because of the pentagon relation, 
trivialize the subbundle satisfying fusion constraints. 

The term 6D is an energy penalty which charges energy 6 for each pair of (5-gon, 
7-gon) pairs created by an F-move (see Figure 2.4a). The number of such defects 
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is counted by the operator D. The expected Levin- Wen 12-body plaquette term is not 
directly included but we will show that it arises at high order by considering the process 
which virtually breaks a r-labeled string, resulting in a pair of "electric" excitation^ 
costing energy 2e, see ifl lll . We consider two triangulations A and A' (and their dual 
nets (or "lattices") and A^') to be equivalent if they are isotopic on 2, i.e. if we 
can slide one to the other. Later, we set S to zero to obtain a lattice liquid, then add a 
string tension term and also tie each net to S in a fixed way (unrestricted isotopy will no 
longer be permitted); but for now we study the crystalline case and need not distinguish 
isotopic nets. The dynamics on the set of nets Nn dual to «-vertex triangulations is 
shown in Fig. 13. II 





Figure 3.1: H-I move, called F-move when coefficients are added as in (1.1). 

Now let us define H^^^ more precisely. First, it enforces fusion rule terms at each 
vertex of each net A^' by penahzing the illegal Fibonacci fusion (see Fig. 13.2b and its 
symmetries. 




Figure 3.2: Illegal Fibonacci fusion. 



Second, it contains terms between states of adjacent nets and A^' which enforce 

, T - Let V, w be the normalized states of H, 



the unitary F-symbol 



shown in Fig. 13. 31 The second terms of //^^ are of the form (id-\v){v\) and (id-\w){w\). 



We discuss the spectrum of //^^ first. //''^ is positive semi-definite and its ground 
state manifold consists of the states ip with {ifr\H^g\t//} = 0. Such a wave function i/' 
is completely determined via the F-symbols by its restriction to a sample net A^*^, e.g. 
a honeycomb. (Importantly, iff is not over determined (frustrated) since the F-symbol 
satisfies the famous pentagon equations.) The ground state manifold may be classified 
according to the number of magnetic particle^ m (of which, in our example system, 
there is only one type). Since we have only imposed fusion and F-moves there is 
no energy penalty for m charges on plaquets, provided that, unlike in the Levin- Wen 



'in the notation of iToll the electric pair may be either (r (gi 1, t Ki 1) or (1 (gi t, 1 (gi t). 
''in notation of flo'l m = r 8i f . 
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Figure 3.3: The states v and iv. Solid lines carry the r particle label and dotted lines the trivial label. 



model, they are allowed to roam ergodically according to the moves (F) which link 
adjacent nets. The magnetic charges on N^^ can return arbitrarily permuted, so the only 
zero energy (unfrustrated) states with y-magnetic charges, j > 2, are the ones that have 
equal amplitude for all positions of the j charges (on all n-vertex nets). In contrast to 
the LW model magnetic charges are localized gapped excitations, the F-moves here 
delocalize them and lead to a gapless continuum of magnetic charges above the ground 
state. In addition, there is also a continuum of gapless "gravity waves," or phase oscil- 



lation across the (not very tightly bound) graph N„ (see section Appendix B i. This is 
analogous to coexisting gapless magnon and phonon excitations in a quantum magnet. 
N„ is regarded as an abstract graph with the triangulations A„ (or nets A^) as vertices 



and edges given by the move shown in Fig. 13.11 The appropriate vertex weighting 



of N„ (see Appendix A i is uniform until the 6D term is added. As explained in 1 
("Introduction") and |3.2| ("Liquid Case") below it is believed that Ai(N„) ~ 

Adding the term 6D induces frustration which causes the ground state wave func- 
tion i/tq to concentrate near the original honeycomb states. When 6 > const ■ y, some 
const * 4 and y the energy scale of the F-symbol, kinetic considerations are over- 
whelmed and there will be a phase transition to exponentially small fluctuations around 
the honeycomb configurations. This concentration alters the weights (see |Appendix A| l 
on the vertex set V{N„) and is expected to gap out the gravity waves. 




Second I - H move here. 



Figure 3.4: One I - H move creates a pair of (5-gon, 7-gon) pairs costing energy S. 

Treating a pair of (5, 7)-gons as the fundamental excitation with cost 6, we see that 
the hexagonal crystal melts (at first order in perturbation theory) for "kinetic energy" y 
associated with the F-move satisfying y > f • 
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Figure 3.5: A second I - H move separates the two pairs. A third I - H move along any of the four bold 
bonds (Figure 3.4) causes a (5, 7) pair to further propagate. Delocalizing the defect at S = Ay in lead order in 
perturbation theory in yjS. 



Now consider a virtual excitation (of energy cost = A) which pulls an electric pair 
(say (r®!, T®!)) out of the vacuum. Because of the nontrivial mutual statistics between 
the magnetic (r® r) and electric (r® 1) excitations, a frustration arises which increases 
the cost of the electric pair if/Y in the presence of j magnetic particles. For small j the 
effect is roughly linear: 

{rf\Hl\>pY) - <«Af I<IC*) - (3-1) 

for some a > and where we have set {^"q" I/Z^^Ii/'q'^ ) = 2e. 

Here a - 7/54 is the energy scale y of the F-symbol constraint divided by the 
number of F-moves required to take one plaquette B of A^o around a neighbor A and 
across an "electric string" (see Fig. 13. 6t . This 1/n scaling of a mirrors that of the 
ground state energy of a one dimensional ferromagnet on a system of length n with 
twisted boundary conditions. 

The splitting which separates the "true" j - Q vacuums now show up at 56th order 
in perturbation theory. (56 - 54 + 2, the 54 counts the steps in Fig. 13.61 to move B 
around A and the 2 comes from first creating then removing the "electric" pair.) The 

true vacuum has its energy lowered schematicall\n by whereas the /' vac- 

uums have an energy reduction of rr. These numbers are each « 55th 

powers of a small and a somewhat smaller number, respectively. One may say that 
perturbation theory predicts, in some regime, a definite splitting off of the true vac- 
uum which, although vanishingly slight, is constant in system size. Thus in summary, 
the timid "crystalline" off lattice approach succeeds in principle but may be judged 
physically useless, because the gap will be tiny. 



^Slightly more accurately by (combinatorial factors) ■ A--y^'*(S + 2e) ^^(2ey 
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Figure 3.6: 9 /-//-moves (or /^-moves at the Hamiltonian level) rotate the three bonds meeting the "dot" 
-1-120°. This rotates hexagon B -60° around its neighbor A from B to B'. 54 moves complete the circuit. 



3.2. Liquid Case 

Let us now move to the other extreme and drop the SD term by setting 5 = 0. Now 
our kets are over a Hquid of "lattices" or "nets" familiar in the quantum gravity litera- 
ture. Let us summarize what is known about the statistics of these nets through theoret- 
ical and numerical study (see Appendix B Appendix D jgl], and |@]). Given uniform 



weight, a weight proportional to total Gauss {- scalar) curvature, a topological quantity 
in dimension 2, or any other known local weight which does not enforce a "crystal," 
the geometry is cactus-hke, with many budding or "baby" universes. These correspond 
to Cheeger constant ^ i2£il (using the combinatorial weights), n - #triangles. That 
is, bottlenecks of size log n are common. In fact, the probability of a bottleneck of 
constant size is » n 2 , i.e. only algebraically small. In a related vein, studying mixing 
times suggests (see [Appendix D| i that A (= Ai of L, see [Appendix A[ ) decays as ap- 
proximately n '^^ when the nets are weighted uniformly. We find a similar exponent 
for the related case of multi-loop rather than net dynamics. 
This seems to present us with two problems: 

1 . Gapless gravity waves 

2. Loss of a length scale 

The first turns out, by itself, not to be a serious problem. It is actually quite inter- 
esting to have a simple mathematical model which manifests gapless modes living side 
by side with protected topological degrees of freedom. In the context of FQH states, if 
the model is taken to be sufficiently comprehensive to include lattice ions, then surely 
their phonons are also an example of this phenomenon. On the other hand, the loss of 
length scale is inherent in declaring kets to be isotopy classes of labeled nets is a seri- 
ous problem. We no longer know if a bond is long or short, straight or wiggly. We view 
with concern the loss of combinatorial protection conveyed by a large regular lattice. 
Recall that on an L x L torus mixing of topological sectors occurs via tunneling along a 



10 



Wilson loop of length L and will be suppressed by a factor of e"<=°°^''^. As noted in the 
introduction, if bottlenecks cause L to be replaced by log L, or even worse a constant, 
then the exponential protection disappears. 

To deal with this problem (2), we introduce a fixed fine-scale lattice on the surface 
S (perhaps writing E as an L x L torus) and regulate our nets to lie within this fine 
lattice. The nets N„ still are restricted to n vertices, n - small constant(L^), but now 
their detailed position in S is pinned as part of the data of a ket |(A, S )); A is regarded 
now as a specifically located, or pinned, n-vertex net in S. We will need to impose 
something that acts like "string tension" that prevents the net bonds from becoming too 
long as measured in the underlying fine L x L grid. This prevents short essential loops 
and so avoids baby universes. As explained below, the bonds become "virtual", only 
their end points are precisely located. String tension can be simulated by establishing 
a hard energy penalty term coB in H, which charges energy co for net bonds longer than 
(max grid bonds (counted by the operator B). Alternatively, a harmonic string tension 
can be imposed. 

Technically the simplest way to incorporate our pinning and string tension terms 
is to alter the basic Hilbert space on which the Hamiltonian is defined. Begin with 
a fine lattice of sites on the surface (such as a torus) and as kets take all pairs: (a 
bond indexed by two sites no more than fn,ax steps apart and thought of as joining the 
sites, a label on the bond). Note that the precise physical placement of the bond is 
not chosen to be part of the data defining a ket0 One may say that the string net is 
"imbeddable" — according to certain rules — but not "imbedded." The bonds at this 
level are "virtual." The labeling just mentioned is from the appropriate set of quantum 
group representations - as is usual - {1,t}. The fusion constraints now specify that 
the virtual bonds first form a trivalent string net and second that the three labels at any 
juncture obey the algebraic fusion rules appropriate to the system of quantum group 
representations being used. An additional "isotopy" terms shifts the location of a vertex 
within the underlying lattice, provided all distance constraints are satisfied. The F 
symbol applies to recoupling virtual bonds. 

Let us explain why equally weighted pinned nets are gapless under these local 
moves. The situation is only a slightly more global version of "the space of all arcs 
transversing a rectangle" -: X. A typical arc will be nearly dense - it will come within 
a constant distance of a positive fraction of lattice points. To define a bottleneck or 
"Cheeger cut" on this space of arcs, consider the mid point m, in terms of arc length, 
of every arc. Let U{L) be the set of arcs for which m lies in the upper(lower) half of 
X. The "cut" is U {^L. Since the probability density if nearly uniform for m in X, the 
Cheeger constant satisfies k ^ j^. Thus by Appendix A | )>= /I. The nets will still be 
gapless after pinning. 

Pinning the net restores the g ^onst i/a scaling for tunneling of quasi-particles and 
hence topological protection. Even if the net is thin (in the y-direction) as in Fig. 13.71 
order L/^max isotopy moves are required to move an excitation around an essential loop 
and so operate on the ground state manifold. 

We will now argue that Hqg - + uB, in the pinned context, supports achiral 



*This simplifies detailed balance for the /-//-moves. 
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4 

excitation path 
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Figure 3.7: A net with a thin part. 



topological phases such as Dfib and the toric code similar to Hiw, the Levin-Wen 
Hamiltonian, yet coexisting with gapless gravity waves. The + superscript in //^^ 
indicates that we have added one more "between fibers" term to our basic Hamiltonian 
H^^, :- "virtual" fusion and "virtual" F-terms. This term raises or lowers the number 
of vertices of the net by 2. The Hilbert space is constrained now to have a maximum 
of n vertices per net rather than exactly n vertices. The new term introduces (removes) 
a "bubble" into a virtual edge: 



with any allowed effect on labels. Rather than write a general formula as before, we 
give as examples the toric code and Dfib cases instead 



toric code: 



V2^ 



(3.2) 



vr 



+ T 



-2 



1 

+ -- 

T 



Dfib: 



Just as we did below Fig. 13.31 the relations of ( I3.2l i are easily converted into pro- 
jector of the form: 

(l-|v)<v|) and (l-|w)<w|), 
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where (in the Dfib case) 



^ unnornmlized 



W unnormalized 




These projectors are the new term in H^^. 

Once a bubble has been introduced in a bond, a succession of 6 (or actually degree 
of face) many F-moves inflates the bubble, carries it around the face and then collapses 
it back to its initial state 




This process is easily seen to enforce the LW plaquette condition: the flux through 
each plaquette is trivial. The question is how strongly it is enforced. Note that there 
is no excitation present in Fig. 13.21 rather we show a circular family of "horizontal" 
terms relating one fiber to another, returning finally to the original fiber Consequently 
the cost of a violation of this emergent "no flux" condition is not a high power of small 
number but rather proportional to the reciprocal of the numb er of horizontal terms in the 
loop. Again think of a one dimensional ferromagnet with twisted boundary conditions. 
The number of horizontal moves is measured by the nets combinatorics, 6 + 2, which 
is 6 for the trip around the hexagon, 2 for creation and annihilation of the bubble. Thus 
despite a suppression by a factor of i, there is a substantial gap to magnetic excitations 
in terms of the bare energy scale of the F-move and the "new" term in H^^. Of course 
the cost of an electric excitation is precisely the bare energy of the "vertical", i.e. within 
fiber, terms which enforce the fusion rules. 

Although in a random net some plaquettes will have more than six sides, the proba- 
bility of s sides decays exponentially with s (see Appendix D i. Thus a small portion C 
of the configuration space N with i-gon plaquettes, s large, has cusp-like geometry (as 
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in the case of hyperbolic geometry) and supports neither small Cheeger cuts nor low 
lying eigen functions. It follows that although a magnetic excitation may be cheap over 
C there is no efficient way to tapper off the amplitude towards zero on N - C where 
magnetic excitations are expensive. The conclusion is that our analysis for s « 6 is, in 
fact, general and shows a gap to magnetic excitations across all the configuration space 
N. 

A detailed comparison of Hgg to the exactly solved Levin-Wen Hamiltonian H^w 
is instructive. The ground states (in the thermodynamic limit) are expected to be bi- 
jective. The excitations of Hgg are, in contrast to Hlw, mobile. To build point-like, 
confined excitations "wave packets" will need to be formed. Combinatorial recoupling 
arguments show that if such packets are confined in potential wells and braided, the 
LW (i.e. Jones) braid representation will be exactly realized (in the strong confinement 
limit). Thus, we may expect that the entire topological structure, the TQFT, represented 
by Hlw is recaptured by H^g. It is true that braiding will excite gapless gravity waves, 
but these are visibly non-interacting with the topological information contained in the 
combinatorics of labeled nets and their recoupling rules. 

We would like to explain more fully this remarkable property of the liquid phase. 
This is the rigidity of topological information maintained in defiance, so to speak, of 
the gapless gravity waves which propagate about. To do this let us speak metaphori- 
cally of the underlying space X„ of (unlabelled) configurations as a "chain". This is 
a reasonable picture since our spectral studies show that the low eigen values of the 
graph Laplacian have inverse power law scaling similar to the l/n^^ scaling of a chain. 
We may very roughly view the quasi-geometry of X„ as a string of length 0{n). We 
should worry that very near the ground state energy we will have states whose topolog- 
ical characteristics "rotate" as we pass from one end of the chain to the other Recall 
our two main exemplars: the toric code and Dfib. Both of these have a 4 dimensional 
ground state Hilbert space (torus) spanned by the states |1), |2), |3), and |4). Imagine 
a system state that is a family of topological ground states that rotates by In as we 
move across the length of the "chain" X„ and so, on the torus triangulation at "chain 
position" X, < X < L, we see the ground state: (cos(27tx/L) |1) + sm(2jTx/L) |2)) . 
Is such a system state a candidate for a low energy excitation as n = L approaches 
infinity? The answer is: "No". To see this look at consecutive "links" in the chain, 
triangulations Ai and A2 with states "Pi on Ai and '¥2 on A2. By the "code property" 
of topological ground states (see fl5l), ^Fi and ^'2 cannot differ by the application of 
a local operator. Passing between ^'i and will cost energy according to the 
term. In fact, this rigidity is quite robust. Up to the usual caveats about perturbations 
inducing exponentially fine energy splittings, it is not possible to deform the ground 
state as one moves through the configuration space X. Since for us the configuration 
X e X is a dynamical variable, this is important. If topological information is stored in 
this novel phase, when it is retrieved X must be sampled. X will be sampled according 
to some distribution and the topological state over the sampled x will then be probed 
by a quantum measurement. The output distribution of our probe , when applied to a 
system ground state, will be independent of the sampled x € X, as desired. 

Because of the pinning and string tension terms, the typical nets in this lattice model 
are qualitatively similar to the boundaries of Voronoy cells produced by Poisson dis- 
tributed centers. We recommend this alternative model to the investigation of interested 
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readers. 

Hqg is not a "lattice Hamiltonian." In particular, it is not defined on a "tensor 
product" Hilbert space (but rather a fiber-wise direct sum of these, one for each net in 
N„). Thus, it is not precise to assert that Hqg is "A:-body" for any k, but it is evidently 
quite simple. One may say that the flux (plaquette) term of Hiw, which is 12-body, or 
more precisely a 6-parameter family of 6-body interactions, has been simulated by the 
F-move, which in these terms is a 4-parameter family of one body interactions. But to 
achieve this, we have resorted to a context where the lattice itself fluctuates and must 
be counted among the dynamic variables. 

3.3. CDT: A(l+l)DHome forAnyons 

Causal dynamical triangulation (CDT) builds layered (l-i-l)-dimensional "space- 
times" by randomly constructing Lorentzian strips as below, where all horizontal bonds 




t pseudo-time = Q 



Figure 3.8: One layer of space-time in CDT. 

have length^ = 1, and all other bonds have length^ = -a < 0, where a is a fixed 
constant, see Fig. 13.81 The action S — /(-^ H- 2A) is the integrated Regge scalar 
curvature (appropriate to Lorentz space) plus a suitable cosmological constant. 

We take pseudo-time 6, periodic. It is known [6] that for suitably chosen a x .7, 
there is a liquid regime of roughly flat Lorentzian geometries on S' x R. This finding 
offers a remarkable solution in cosmology to the persistent problem of baby universes. 
We recommend for further study the possibility of importing this innovation into con- 
densed matter physics. The same action can be used to define a density on string nets 
supported neat flat geometry, and so preserve topological protection in anyonic models 
based on these geometries. 

To visualize the braiding of anyons, described next, picture S' x R as - (0) via 
(6,p) i-> (eP, 0). The geometry of R^ - {0), i.e. its Lorentzian triangulation, is explicitly 
among the dynamic variable, but in addition the bonds of the triangulation are labeled 
from a (quantum group) label set, which in this paper is {1,t). "Singularities" of the 
labeling (as explained in detail in Fidkowski et ai, see [10] ) - annular regions where 
the state cannot be extended over the disk to a vacuum state - are the "quasi-particles", 
or anyons, of Dfib. So a loop of states is a loop of annular Lorentz geometries together 
with anyons. 

One might wonder when importing a 2-dimensional net model from a ( 1 + 1 )-dimen- 
sional quantum gravity model, whether "causality" in the model will prevent braiding. 
If information is not allowed to flow backwards in "time" (pseudo-time = 6) we might 
be unable to braid anyons since they can only move forward in the 6 direction. This 
may appear to limit their possible braidings, but in actuality it does not. A full counter- 
clockwise 27r-turn generates the center C of each braid group, B„,n > 3. Thus as the 
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anyons move radially back and forth, their overall progression in the pseudo-time = 6 
direction only multiplies the braid by a central element - corresponding in each irre- 
ducible sector of the Jones representation to an irrelevant overall phase. Consequently, 
the causality of the construction (in pseudo-time - 9) does not restrict the image of the 
braid representations and building topological phases is not hampered by a causality 
constraint. 




«— f = time 



6 = pseudotime 



Figure 3.9: Anyon trajectories superimposed on a central twist. 



4. Conclusions 

Because nuclei are heavy lattices in condensed matter are generally thought of as 
fixed or classical degrees of freedom. It is true that in chemistry superpositions of iso- 
mers can be important but generally the lattice is not takes to be a dynamical quantum 
mechanical variable. In this paper this is exactly what we do. We have not forgotten 
that nuclei are heavy; we imagine that there may be models in which some electron 
degrees of freedom define a lattice and others decorate it and that both should be al- 
lowed to fluctuate. This paper is not about a specific model of this kind but rather a 
preliminary survey of the hazards and prospects that await us "off' lattice". Our focus 
has been entirely on building topological phases, although off lattice models may have 
wider applicability. 

We have seen that the chief hazard is uncontrolled fluctuations in the now random 
geometry of the lattice, a phenomenon colorfully called "baby universe" in the quantum 
gravity literature. These fluctuations threaten to destroy the ratio of scales L/a, system 
length / lattice spacing, on which topological protection, error ~ exp(-const.L/a) de- 
pends. We have also seen our three approaches (methods of birth control) to retaining 
topological protection. Briefly they were: 1) a minimally fluctuating crystalline phase, 
unsatisfactory due to a vanishingly small excitation gap A for the topological phase. 2) 
Pinning the fluctuating lattice to a background (lattice or continuum). This seems to 
work but sacrifices some of the simplicity we hoped to find in off lattice models. 3) 
Causal Dynamical Triangulations (CDT). Here we borrow the solution (as well as the 
problem) from the quantum gravity community. 



In order to evaluate the impact of geometry fluctuations, analytical ( Appendix C i 



and numerical (Appendix D i work was done on the statistics of loop gases and string 
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nets. In Appendix A and Appendix B Cheeger's theory relating the geometry of man- 
ifolds to their vibrational modes is adapted to the infinite dimensional context to con- 
struct estimates of the spectrum of Hamiltonians derived from our models. In particular 
we develop a method for constructing upper bounds to the spectral gap of a Hamilto- 
nian H by Monte Carlo studies (of both the gap and estimates for Cheeger's constant) 
on a weighted graph G derived from H. We find (method 2 and appendices) regimes 
in which topological information is protected while coexisting with gapless "vibra- 
tional" modes across the space of geometries. Such results encourage us to regard the 
off-lattice approach as viable and worthy of continued investigation. 
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Appendix A. Spectrum of Graph Laplacians and Other Local Hamiltonians 

We start with the following data: a finite dimensional Hilbert space 'H spanned by 
a preferred set of basis kets a Hamiltonian H . "H, and a known ground 

state wave function i/^o = 2,- fl; |/) for H. We construct from the data a weighted graph 
G whose Laplacian X. . 'H ^ 1-( is easy to study numerically. We focus on the first 
eigenvalue /li(-0 and verify the Cheeger inequalities for G: 

2hc>Ai(£)>^. (A.l) 

If it should happen that Ai{X.) — » (as a scaling limit is taken), we may also conclude 
that the original H is gapless (in the same limit.) This is because a small Iiq means a 
neck in the set of kets {|/)) with little coupling from cH (where c is a positive running 
factor, perhaps proportional to (system s/ze) arising in the proof) from left to right 
sides of the neck. The trial wave of the form if/i = /'ii/'q* - ^2i/'o^'''(^b ^2 > 0) will be 
orthogonal to i/tq and satisfy c({i//i \H\tfr[}- H li/ro)) 0. If this rate of convergence 
to zeros is faster than c, it will imply H gapless. Conversely, if we know a quantum 
mechanical system is gapped (e.g. the Levin-Wen model [3]), it will imply a tightly 
connected geometry for the appropriate weighted graph of string net configurations. 

The reader may wonder what good is a method for studyinig the "gap" if it requires 
knowledge of the ground state (/tq. In the case of a topological phase, one may begin 
with a formula for the ground state wave function (given by d-isotopy, or a chromatic 
evaluation) and from this, attempt to build a gapped Hamiltonian. This appendix pro- 
vides ammunition for shooting down such Hamiltonians (i.e. showing them gapless) 
as in iflil . 

Here is the construction. The vertices of G are simply the index set {/) for the kets 
of H. We set the edge weight w,j = c|//,j| provided / 4^ j, for a positive constant c yet to 
be determined. The wave function weight at / is dj = |a/p. Write di - cY^ji^i \Hij\ + w,-,- 
where c is the largest (/ independent) constant allowing all w„ > 0. This fixes c and 
the w„. The cw„ are to be thought of as weights on loops at /. Thus, G has edge 
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weights cwjj and vertex weights t/,. Following Chung fs*], there are "unweighted" and 
"weighted" operators L and X., both of which are symmetric and have a zero mode: 

\-Wij otherwise 
£ = T-^^^LT-^^^, (A.3) 

where T is the diagonal matrix with T,, = c/,-. Explicitly, 

fl-^ if/ = ; 

= \ -'V,/ ' ,u • (A.4) 
_— otherwise . 

Example A.l. As a sanity check on the method, we check explicitly, in a toy model, 
that the choice of basis only affects the spectrum slightly (by a factor of 2.) We explore 
in the simplest case the dependence of the spectrum of £. on the choice of bases for TY. 
Let-K^C^and 

/ cos 9 - sin 6\ 
l-sin0 -COS0/ 



e 



Hg annihilates (^^'^^^21 ' need to solve, with c as large as possible, the equations 

, e 

di - sin - = c| sin0| + wii 
d2 - cos 2 ~ ''I ^™ ^1 

Recalling that sin0 - 2 cos | sin |, the natural (and correct) guess in the interval < 
< ;r/2 is c = i tan f . This yields: 

wii = and W22 = cosO 

Substituting, we find: 



1 tan I 
tan I tan^ | 



Solving for the eigenvalues we obtain Ai = 1/ cos^ |, which in the considered interval 
< < 7t/2 varies only between 1 and 2. In the other intervals we get similar results. 

With this small check of quasi-invariance of Ai under basis change, we derive the 
Cheeger inequalities in the relevant weighted graph setting, closely following [SJ. 

X, acts on functions of G by (left) multiplication. The lowest eigenvalue is /lo = 
with eigenfunction fo(i) - -\/di. When G is connected, /Iq is non-degenerate. We will 

if-Cf) 

be concerned with the next eigenvalue Ai = inf — computed with respect to fj., 

the measure or vertex weighting with weight(/) = dj, for / orthogonal to constants. We 
denote Ai by Ac or just A and use / for its eigenfunction. 
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We define the Cheeger constant h: 

h - min _ , 

s^v rmn{Vol{S),Vol{S)) 

where 5 is an arbitrary subset of the vertex set V{G), S isV\S,F„ denotes the weight 
of edges between S mdJ,F^= ^ . Finally, VoliS ) = ^ di, VoliJ) = ^ dj. 

ieSJeS '^^ jeS 

Theorem A.2. 2h>A 

Proof. Let S achieve h and set a = VoliS ) and b = Vol{S ). Define a "trial" eigenfunc- 
tion: 



— la' 



fi 

We have, from the Rayleigh-Dirichlet integral: 

A<F{S,S) " ' 



FiS,S) 



1 M 



[a 

2FiS,S) 



mm(Vol{S),Vol{S)) 
2h 



□ 



Theorem A.3. /I > f . 



Proof. Let functions f,k : ViG) -> M be related by A:, = dif. Now the Rayleigh 
quotient: 

becomes A when minimized among fc, orthogonal to c/,, equivalently by /, orthogonal 
to constants. We assume / is such a minimum. Thus XJc = Ak. 

We index the vertices i of G in /-increasing order, f < f+i, and without loss of 
generality assume ^ - ^ '^j- For i e V let c; = ^ wj^, measures the j'th 

fi<0 fj>0 i<i<k 

— Ci 

"cut" between 5, - {j < i] and S . Set B - min z . Clearly B> k. 

'^y mmiVol{S),Vol{{S))) 
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We set = {i\fi > 0) and the set of edges with at least one endpoint in V+. 

„ f fi iff ' £ V'+ - , 

Finally, set g, = ^ q otherwise ' compute A: 



y fi {Zii,j)eE, »'ij(fi - fj)) 

(cutting off some numerator terms) 
Z Wjjigi - gjf 

(Z(f.J)E£ Wjjigi - gjf) {lj(i,j)eE Wjjigi + gjf) 
{Zievgjwi) {ZdjyeE Wijigi + gjf) 

(by Cauchy-Schwartz) 

{Z(iJ)eE Wu(gf-gj)f 
{LieV gjw,) {ZdjieE Wijigi + gjf) 

(Since Mgi + gj)^ ^ 2 2 g]wi) 

E V 

(discarding cross terms from the numeral 

{i:iCi\g] - g]\f {i:iPyoi(Si)\g]-g]\f 



(telescoping the sum) 
_ {T.iP^igff _/5\h^ 
^2{Zygf4^ ' ^ 2 

□ 

To better understand the proof of Theorem lA. 31 we summarize Cheeger's original 
argument in the context of a Riemannian manifold M. Define Cheeger's constant h by: 

, . _ Area(S) 

h = mt 



S separating M Volume(M) 
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Let / be the first eigenfunction of the Laplacian A orthogonal to constants. 

JfAl Imlf\ {f\fm\f 
~ If ~ If If - {fff 



4 

Define t - as a parameter on M and apply the co-area formula to the f-levels to 
obtain: 

J (V/) = J Area(t-\eve\)dt <h J Vol[0, t]dt 

t-£-dt ^-h I tdVol = -h I fdVol. 

Thus, 

^ 1 {-h JfdVolf _ h2 

□ 

We offer a protocol which may succeed in verifying that a quantum mechanical 
Hamiltonian H„ : 'H„ — > 'H„ is gappless above its (known) ground state (i^o),, as a limit 
n — > oo is taken. 

Protocol. Select preferred kets |/) for (we do not clutter the notation by showing the 
dependence of the index set {/) on n.) use H„, (iAo)n to construct the weighted graph G„ 
as above, and set /I" - Ai(£.c,}- Recall that the construction of G„ requires extracting 
a constant c„ (in our two dimensional example c" - - (2 cos |)"',) the minimal 
suppression factor for interactions Hjj requked to normalize the vertex weights di :- 

|a,p be positive and with 2, di - 1. Compute the ratio ^"f^ , where D„ - max, |a,;„p 
for i/fQ.n = Zi ai,„ \i). 

Claim. // '^"^ — » 0, then H„ is gapless, i.e. Ai(H„) — A(){H„) — > a.? n — > oo. 

Proof. We have checked A" = /1i(-Cg„) ^ -f^, so /z„ = Iiq^^ < ylA". As in the proof of 
(IA.2l l. let S" achieve h„ and define: 

r^i)) = -(- + ^\ i e S'\a„ = Vol{S") 

rim = -U- + i 6 5^", bn = Vol(S"). 
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Then, suppressing the n super/subscripts, 



<<Ai|//|<Ai)-<'Aol//|'Ao> = £1 -£() 

< y--^\ai\\aj\\Hij\ 

ieS 

< --^DFJS,S)c-' 
a I - a 

=: DFc-\ 

a 1 - a 

where a - Vol(S) and without loss of generaUty, Vol{S) < Vol{S). From Cheeger's 
inequahty iA.3\ . we have: 

F I— 
- = /! < VIA, 
a 

so 

El - Eq < -DFc-^ < D V2lc 
a 

□ 

This protocol allows a systematic approach for vetting models which produce known 
topological wave functions as the ground state (say on a 2-sphere) but may not be 
gapped above the ground state. There have been previous successes in showing models 
gapless by finding directly the Cheeger cut into V(G) - S US 1 14]. The present proto- 
col may be more practical as less geometric insight is required, unfortunately plugging 
in the analytical bounds from Appendix C into the claim we find ^ ^ n'l\ Using 

the numerical scalings in Appendix D we find ^" ^ x which does not approach 
zero as n — > 00 either Further geometric insight into the graph G„ might allow one 
to use A" rather than yfJ" in the claim, yielding n"^^"*. This would be legitimate if G„ 
looked spectrally more like a tree than a line. 

The d = 1 loop gas 1 14] has the surprising feature that the very same ground state 
arises as a gapped and gapless ground states of two different Hamiltonians |il3j] . A 



second example (d — V2) was proven OlSll via decay of spatial correlators never to 



arise as a gapped ground state for any local Hamiltonian. However, in many cases 
unlike the above, one will not be so fortunate to find a narrow cut for G. Rather, more 
generically one may expect to learn something about the spectrum (A) of X on G and 
perhaps some properties of the first eigenfunction / via Monte Carlo methods applied 
to G (since this problem is completely classical). In this case, one should try to use the 
protocol. In a gapless system, to find the precise power at which A" — > 0, more refined 
trial wave functions involving a gradual, not abrupt, phase change across the cut should 
be studied, as in [ 14]]. 

Note A. 4. Our protocol can be used in contrapositive form to argue that the A associ- 
ated to certain weighted graphs of configurations cannot decay too quickly in system 
size when we know that the (weighted) graph arises as the ground state of a gapped 
Hamiltonian, such as the Levin- Wen model. Specifically, when applied to the Fi- 
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bonacci anyons, one may argue that the set of subgraphs G of the honeycomb when 
weighted by its topological evaluation (see 1 10]) and supplied with edges correspond- 
ing to a bounded number of F-moves and local circle creation/deletions on a finer scale 
honeycomb must have Aq decaying no faster than (system size)"' . 

Monte Carlo methods may eventually be able to extract some information of the 
eigenfunction / associated to /Ii(Xg")- It is reasonable to suppose that knowledge of 



/ could refine the previous protocol. The final paragraphs of Appendix A present, 
schematically, the outlines of a complementary approach to extracting information on 
the quantum mechanical spectrum (//) from the first eigenfunction of the classical £.g- 
Let us use / to build a trial wave function *Fi = /^Fq from the ground state *Po of H. 

M,H - M).H < <^i|i/|<Ai) - (MHlifro) = ^<.M)|r„|.M)) - ^<tAoir„|.Ao) (A.5) 

a a 

where we have written H - 2„ Tq, as a sum of local terms. For any term Ta which 
acts at state / e V(G) we should study the variation of the quadratic forms on the 
right hands side of equation lA.5 1 at second order in the gradient V/. (The 0-th order 
variation vanishes since / is normalized, J f^dfi - J l^dfi - 0, Xii ldiP = 1- After 
summing over a, first order variation must also vanish since {i//o\H\il/o) is critical for 
(actually minimizes) expectation.) 

At 2nd order in V/ and with a fixed, 

r.h.s.„(lA3 z A^^'WVfif (A.6) 

where /l™^" measures the largest eigenvalue of Ta after normalizing all eigenvalues to 
be positive. In ( IA.6I ). / ranges over states on which r„ operates. If n is the maximum 
number of terms Ta operating on any state |/), we may "integrate over a" to obtain from 

(KB-- 



rh.s.Q < n^r" ^ ||V/;||2 

ieV(C) 

< riA^'^^'Ac, for ^f'' = max{^™ 

So, at least schematically, there should be an estimate: 



We finally wish to mention a related paper 01611 which uses similar methods to argue 
for the existence of gapped models. 



Appendix B. Outer Planar Triangulation 

An outer planar triangulation (OPT) is a triangulation of the n-gon P„ in which no 
new vertices in the interior disk are permitted. The «-gon is given a fixed base point 
vertex and orientation. Thus, for n = 3, 4, 5, 6,. . . the number of OPT are 1, 2, 5, 
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14, In general, \0PT„+2\ - Cn - ;7fr(^"). the nth Catalan number. This statement 

may be familiar as the correct counting of dual planar trivalent trees. 

Let G„ be the abstract graph with vertices OPTn and edges determined by "di- 
agonal flips" defined on quadrilaterals made from a pair of triangle sharing a bond. 
Give all vertices and edges of G„+2 unit weight. The spectrum L of G may be simi- 
lar to that of the more interesting case of triangulations of the 2-sphere and because 
there are simple asymptotic formulas for the Catalan number, we can explicitly com- 
pute a lower bound kn < 0{n'^^^) for our Cheeger-like isoperimetric constant k„ = 

min — ^ — ^— where E(S,S) counts edges from 5 to 5 = V(G„) \ S 

scV(c„)inin(yo/(5),yo/(5)) 

and VoliS) = # vertices in S . We thank Oded Schramm for guiding us through this 
example. 

For simplicity (only) take n odd. Now there will certainly be a unique "central 
triangle" A with the property the three connected bits of sides(P„) all contain less than 
I sides. Call the "lengths" of these three bits j > n\ > n2 > ns, n\ +n2+n3 = n. We 
divide OPT„ into two disjoint pieces, thick„ U thin„ - OPT„, and thick„ n thin„ = 
according to whether > ^ (called thickn) or «3 < (called thiun). 

We will use the well known relation Cn±const 

= 0(l)4"n 2 and in the future use » to 

absorb the 0(1). 



\thick„\^ 4'"n^h''^n^H"'nl 

§>ni>n2>«3aTO 

4"(nin2«3)"^ 

§>«i>n2>«3aTC 

^ 4"(n3)-i 

0(«2) terms 

« 4"n-i 



and 



\thinn\~ ^ (#(ni,n2) with n3 = i)(#(ni,n2, 13 = i)) configurations 

i<.v<fji 

~ *(4"'nj^^4""n2M"^i"5j, for typical ni,n2 with ni -I- n2 
= 4"n-'^ 

So both thick and thin portions of OPT„ have 0(1) proportion of all the vertices on 
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G„. 

Now consider the probability of being "near," within a of the boundary between 
thick and thin: Irij, - < a. Call such configurations "boundary" or da, \da\ ~ aA"n^'- 

_ 9 

as there are Oina) such numerical configurations each occurring order 4"n 2 ways. 
Only diagonal flips on one of the three sides of the central triangle can possibly aff'ect 
membership in thick„ and thm„, and we should estimate how many such flips can relate 
thin to thick. The largest contribution comes from flips on side n\ (or equivalently nj) 
in which a vertices of the «-gon move to «3 where a - o(n). We estimate the number 
of such G-edges as follows: 




Figure B.l: 



Ei{thin„, thickn) ~ constn ^ (4"n 2) 

Above, const n reflects a summation over a. Configurations with n^, - jq - a are 
counted in the first term; the fraction crossing from thin to thick upon the flip indicated 
in Fig. IB. H is given by the second term within the sum. 



b - b) 



3 



E\{thin„,thickn) ~ 4"constn 



4" const n(n -n 2) 



Neglected terms, such as a comparable to n, are down by a power iT'- and have 
been dropped. 

Putting the three calculations together, we conclude that the isoperimetric Cheeger 
constant k satisfies 

4"n-3 



k ^ 



4"n"2 

This means that the valence normalized Cheeger constant h appropriate to random 
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walks list] satisfies 

_ 3 

h ^ n ^- 

and that the mixing time is at least Oirr'^^). A mixing time of a: n corresponds to the 
usual graph theoretic notion of an "expander." Numerical investigation of this model 
indicates that the actual mixing time is 0{n^). 

Appendix C. Analytical results for the off-lattice loop gas 

Appendix C.l. An off-lattice loop gas model 

In this appendix, we study analytical properties of an off-lattice loop gas model. 
The basis states of the model are configurations of non-intersecting, indistinguishable 
loops, identifying loop configurations related by isotopy. Loop configurations with 
at most loops can be represented by unlabeled rooted trees with at most A' nodes, 
excluding the root node. Using the recursion relations of Ref. {vf\ the number of such 
trees (i.e., the number of loop configurations) for a fixed number of nodes n is given by 

C(n) = y C(n - /t) y mC(m) , (C. 1) 

n — 1 ^—^ ^—^ 

k=\ m\k 

where n > 1, C(l) - 1, and "m|A;" denotes all m which are factors of k. A similar 
expression exists for the number of leaves (excluding the root) of unlabeled rooted 
trees, 

«-i 

L(n)^Y,C(n-k)J]L(m) (C.2) 

k-l m\k 

withL(l) = 1. 

The Hamiltonian of our off-lattice model acts locally by the three types of moves 
shown in Fig. IC.ll a) The inflation move corresponds to creating or annihilating a 
loop, b) The surgery move is merging of two loops, c) The self-surgery move is a 
surgery move of a loop with itself. We define the Hamiltonian // as a sum of projectors 
performing inflation, surgery, and self-surgery moves such that for the ground state 
wave function li^o) we have //|i/'o) - 0. The ground state wave function then becomes 
an equal-weight superposition of all loop configurations / 

m^Yj^l). (C.3) 

Note that the Hamiltonian takes the form (up to rescaling) of an unweighted graph 
Laplacian L (see Appendix A): each transition (via inflation, surgery, of self-surgery 
moves) from a state |a) to a state \/3) gives an entry of -1 in the Hamiltonian matrix, 
and the diagonal elements are H„a - - Yipi^a Haj3, is., the diagonal element count 
the number of transitions out of state |a). 

Topological protection 

We consider the loop gas on an annulus (periodic boundaries in one direction). A 
particular loop gas configuration can be represented by a tree where one leaf marks the 
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Figure C.l: Loop gas moves. From top to bottom: a) inflation move, b) surgery move, and c) self-surgery 
move. 



inner edge of the annulus and the root corresponds to the outer edge of the annulus as 
illustrated in Fig. lC.2l The surgery move alters the number of loops that wind around 
the system by +2; the parity of the winding is hence a conserved quantity. 

Is it possible to locally distinguish even and odd winding sectors? In a lattice 
realization of a loop gas, such as the toric code fP], the expectation values of any local 
operator in these sectors split by at most an exponentially small amount - the hallmark 
of topological protection. In an off-lattice model, on the other hand, the splitting of 
these winding sectors turns out to be only algebraically small. To see this, consider the 
average number of leaves in a sector with parity p. The difference between the odd 
and even winding sectors 



can be computed using Eqs. iC.2\ and (IC.lb and is found to decay algebraically as 1 /N, 
which is also illustrated in Fig. IC.4I 

Appendix C.2. Proof of gaplessness of loop gas Hamiltonian 

A loop gas configuration of n loops can be represented by a rooted unlabeled ran- 
dom tree with n + 1 nodes. In the following, we shall refer to this tree representation. 
The ground state of the Hamiltonian is given by 



(C.4) 




(C.5) 



27 




Figure C.2: Loop gas on an annulus. 




Figure C.3: Height-changing moves (inflation and self-surgery — surgery moves do not change the height) 
are only possible if a tree has only one leaf at its maximal height level, like the configuration shown here. 
Such configurations then allow for one inflation move and one self-surgery move that changes the height 
by 1 (left panel). In contrast, there are many self-surgery moves possible that change the height by 2 (right 
panel). 



where the sums runs over all possible rooted unlabeled random tree configurations, and 
Cn is the total number of such configurations with at most nodes, 



(C.6) 



n=l h=0 



where C(n, h) denotes the number of tree configurations with n nodes and height h 
where n> h. The action of the Hamiltonian is such that 



H\il,o) = 0. (C.7) 

The terms of the loop gas Hamiltonian (inflation, surgery, self-surgery) are equivalent 
to the following modifications in the tree representation (we hst only the direction of 
the moves that remove a node): inflation corresponds to removing a leaf of the tree, 
surgery corresponds to "fusing" two sibling nodes into a single node, and self-surgery 
corresponds to "fusing" a child node with its grandparent node (all children nodes 



28 



0.016 




0.01 0.02 0.03 0.04 0.05 

1 IN 



Figure C.4: Power-law decay of A(N), the difference of the average number of leaves of unlabeled rooted 
trees in even and odd winding sectors for an annulus. 

of the child node become sibUng nodes of its parent node, the whole process can be 
visualized as "folding" the tree such that child node and grandparent node become 
one node). The different moves in the Hamiltonian change the tree height by either 
(all surgery moves, some self-surgery moves, some inflation moves), 1 (some inflation 
moves, some self-surgery moves), or 2 (some self-surgery moves), and thus 

-{a\H\a) = ^ (J3\H\a)+ + Yj ^^"^"^ 

We define a "constrained" number of configurations at height h, 

C(h):^ Y C{n,h), (C.9) 

n, where g(n)>h 

where g(n) is some function of n to be defined below. Using this definition, we make 
sure that only a constrained number of configurations is included. Next, we consider 
the state 

h=h ^ ' {a, h„=h, g{,n„)>h\ 

where m > \,h - Int[fe(A^)] (Int[ji[:] denotes the smallest integer number larger than 
x), and C - 1/C(h). The following proof relies on the inequality C(n, h + I) < 

C{n, h) that should be valid for all configurations included in the above trial state. This 
inequality can be satisfied by proper choice of the functions g{n) and k{n), see below. 
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The states li^o) and are two distinct orthonormal basis states, 



-j^ mh-l -j^ -j^ mh-\ ^ 

h=h ^ ^ " {a, K=K g{n„)>h} h=h ' ' 



and 



mh- 1 

z 



expilnih / /;) 



-j^ mil- 1 

— ^ exp(2mh/h) - 



2 <«i^> 

{(y, h„=h, g(nct)>h} 



cc 



h=h 



The energy gap can be estimated as 
A£ < <iAi|//|<Ai> 

-, mh-\ mh- 

= ^ZZ 



1 "v-i' "v^' exp(27r/(/z - h')/h) 



C(h)C(h') 



h=h h'=h 

Using relation (IC.8i t, and that {p\H\a) < if a ^ ^, we obtain 

m/7- 1 



^ ^ </3|//|a). 

(j3, /!^=/i', g(n,,)>/i'| Iff, hc=h, gM>h] 



AE < 



^Z 



h=h 
mh-l 



2 cos(27r//i) 



jZ 

h=h 



(C(/z))2 C(/i)C(/z+l) 
2 2cos(27r//i) 



(C(/!))2 C(/i)C(/z + 2) 



2 J] |</3|//|ff)| 

1/3, his=h+l, g{nii)>h+l] {a, h„=h, n„=«^-l) 

^ 2 l<;0|//|a>|. 

1/3, hi)=h+2. g{nfi)>h+2] {a, /i„=/2, «„=«/j-l) 



A state I/?) of height /z has at most one inflation move transition to only one of ah 
states \a) of height h - I. The same appHes to self-surgery transitions that change the 
height by one. A tree of height h can have at most h - 2 self-surgery transitions that 
decrease the height by two. Using these estimates, and that C(h -H 1) < C(h) (since 
C(n, h+ I) < C(n, h) for all configurations included in the trial state), we obtain 



1 

AE < —a\ sm 
C 



mh—l 



h=h 



2C(h -Hi) (h- 2)C(h + 2) 



(C(/i))2 



(C(/z))2 



, mh- 1 , 

— a2 sin (n/h)hm ^ ~ ^2Sm(7T/h)mh — — 

h-h 



for large h, where a\, 02, and 03 are constants. 

We need to choose the functions g(n) and k(n). The authors of Ref. ifisll showed 
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Figure C.5: Growth of H' (value of height where C(n,K) is maximal at a given n) as a function of n as 
obtained from l IC.lU . The growth is proportional to V"- 



that for n — > oo, 

C{n,h) ~ Cin)2bJ—0* V k\2k^/3^7T^ - 3) exp(-A:V/32), (C.ll) 
y n j—f 

k>l 

where p ^ 0.3383219, b ^ 2.6811266,/? = 2 ^fUlhb^Jp, andlfllllTl] 



2 V''!' 

Eq. (IC.l lb is asymptotically valid for arbitrary but fixed 6 and {6 ^J\ogny^ < h/ ^/n < 
6 -y/log n. It can be seen from Fig. IC.5l that the number of configurations at a given n 
is largest for h'(n) - 2.1 yfii. It follows from Eq. dC.llb that C(n,h + 1) < C(n,h) 
in its region of validity if h > h*{n). It is easy to check that C{n,h + 1) < C(n,h) 
for any configuration in Eq. (IC.lOb if we choose g{n) - 6 -y/nlogT? and k{n) = 2.1 ^fn. 
Here, 5 is chosen in such a way that there exists, for a given A^, at least one n such that 
h < h < g(n). Using this choice of g{n) and k(n), we obtain the following estimate for 
the energy gap 

AE < (C.12) 

where as is a constant. 

It is likely that C{n,h + 1) < C(n,h) as long as h > h*(n) (and not only for 
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6 ^jn log n > h> h*(n) as in the previous paragraph). If this is true then we can derive 
a tighter bound for the gap. Indeed, asymptotically h > h*(n) for any h in Eq. dC.lOb if 
we choose the following function k{n) - an", where 1 /2 < /c < 1 and a is a constant. 
In this case, the gap scales as A^"*^ and the upper bound is A^"' . The states with k that is 
close to 1 have support only on an exponentially small number of tree configurations 
and thus detecting the A^"' scaling in Monte Carlo simulations seems unfeasible, see 



Appendix D 



RescaiingJ this Hamiltonian to a graph Laplacian X , we obtain a scaling of the 
gap as A^"^^^ (from Eq. ( |C.12t ). Plugging this into the gap estimates for general local 
Hamiltonians in Appendix A we find an upper bound of VA'^- We will see below, in 



the numerical results of Appendix D that the gap actually scales as A^ , which is 



still not enough by itself to prove gaplessness of any model. 
Estimating Cheeger's constant 




Figure C.6: Scaling of the average number of moves that take a tree across the mean height ft*. The perfect 
linear behavior shows that the number of moves indeed satisfies the linear upper bound. 

Since the gap of the graph Laplacian seems to be too weak a bound for the Cheeger 
constant, we next try to estimate the Cheeger constant directly. Using the same cut as 
in the above proof, we consider a split of the configuration space into those trees which 
are smaller or larger than the mean height h* . The fraction of trees at this boundary 
can be obtained from equation ( IC.llb to be C(n, h*)/C{n) ~ 1/ V"- Multiplying this 
with the number of moves across the cut, which can be bounded by 0{h*) ~ 0( yfn) we 
obatin as estimate for Cheeger's constant 

h < 0(l/^/N)0(h*(N)) < 0(1/ Va^)0( Va^) = const., (C.13) 
which is now border line regarding the absence of a gap. We hence tried to check 



'We drop an exponentially small fraction of states, corresponding to some trees whose height scales 
slower than 0( V") ™d whose connectivity scales faster than N 
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numerically whether the number of height changing moves, which we have bounded 
by h* might grow slower then linear However it turns out, as is shown in Fig. IC.6l that 
the scaling indeed satisfies this bound. 

Appendix D. Numerical study of ofF-lattice loop gases and string nets 



Figure D.l: String net moves. From top to bottom: inflation move, /-move, and surgery move. 

In this appendix, we turn to a numerical analysis of the ofF-lattice loop gas and 
string net models by Monte Carlo and exact diagonalization. The off-lattice loop gas 
model has been introduced in the previous appendix. In a similar fashion, off'-lattice 
string nets can be defined as indistinguishable (unlabeled) planar trivalent graphs, where 
we exclude configurations with bubbles or parallel edges. We define the system size 
as the maximum number of faces, which is related to the number of vertices via 
- (n,, + 4)/2. The Hamiltonian again takes the form of a graph Laplacian and is 
defined by the three types of moves illustrated in Fig. ID. II The string net ground state 
again is an equal-weight superposition of all string net configurations s 



The definition of the graph Laplacian X. in Appendix A includes the constant 
c - min,(ii,7 Y^ji^i \Hij\). The sum in the denominator is basically the (weighted) ver- 
tex degree of vertex (basis state) /. It typically grows faster than the system size A^. 
However, the number of vertices (or basis states in the Hilbert space) for which the 
weighted vertex degree is not bounded by a linear function of is exponentially small. 
As an example, surgery moves within a plaquette of a string net grow like the square 
of the the number of edges in the plaquette, since one can do surgery between any pair 
of edges. Since plaquettes with a large number of edges are exponentially suppressed 
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Figure D.2: Distiibution of the plaquette sizes in the off-lattice string net 

(see Fig. ID. 21 ). we discard those exponentially rare states and restrict the Hilbert space 
to 'H' such that a basis state |/'> belongs to iff |/') € and Y^ji^i' \Hi'j\/dr is 
bounded by a linear function of A'. The graph Laplacian X.' is then defined following 
[Appendix A[ with c' proportional to ^ 

The graph Laplacian is gapless by definition. In the following, we also demon- 
strate the gaplessness of NX,' for both the off-lattice loop gas and string net models by 
numerically determining the gap to the first excited mode of X' ■ 



Monte Carlo method 

One can extract the gap of the graph Laplacian from classical Monte Carlo simu- 
lations 1I2O] by ensuring that the Monte Carlo transition matrix is proportional to the 
graph Laplacian X''- 

T ^(\-a)l + aX', (D.l) 

where T is the transition matrix and a is the coefficient of proportionality. 

We perform Monte Carlo simulations by first calculating the number A^moves of pos- 
sible moves for a given configuration. Let A^moves - l/c'. Then A^moves is a linear 
function of A^ and it is larger or equal to A^moves for any configuration. We randomly 
pick one of the possible moves and accept it with probability A^moves/Mnoves- The total 
probability to make a move is 1 /Nmmes and it is the same for any move. The probability 
to stay in the given configuration is 1 - A^moves/A^moves- This transition matrix is equal 
to the graph Laplacian X' (a = 1). 

The enumeration of possible moves can be implemented very efficiently for rooted 
trees allowing us to access large system sizes. However, this is not the case for the 
string nets as one needs to check for graph isomorphisms for every possible move, 
which restricts us to considerably smaller sizes. We use the isomorphism test suggested 
inRef. jllll . 

The gap is related to the autocorrelation time ta of some observable A as 
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Figure D.3: Exact diagonalization results: The gap of the graph Laplacian as a function of the inverse 
system size i/N obtained from exact diagonalization. Results for the off-lattice loop gas are shown on the 
left (1/c' = 3200^), and for the off-lattice string net on the right (1/c' = 30^). 



where ta is measured in Monte Carlo time. The observable A must be chosen carefully 
- it must couple to the lowest mode in order to extract the gap. 

Exact diagonalization results 

We first analyze the spectrum of the graph Laplacian for off-lattice loop gas and 
string net model using exact diagonalization. In particular, we calculate the lowest gap 
using the Lanczos algorithm [22] as shown in Fig. ID. 31 For the off-lattice loop gas, we 
find that if we only consider inflation and surgery moves, the graph Laplacian times 



the system size is clearly gapless - consistent with the proof in Appendix C Adding 
self-surgery moves the gap of NX.' appears to extrapolate to a finite value. The same is 
seen for the off-lattice string net. However, this apparent convergence is misleading as 
we will see below in Monte Carlo simulations of larger systems. 

Monte Carlo results 

To determine the gap of the graph Laplacian in Monte Carlo simulations we mea- 
sure the autocorrelation function of the tree height (for the loop gas) or graph diameter 
(for the string net). As shown in Fig. ID.4I we find that, for small system sizes, the 
autocorrelation functions couple to high energy modes resulting in a fast initial decay 
before turning to a slower asymptotic behavior corresponding to the smallest gap. This 
makes it difficult to extract the gap for large system sizes, since at long times the au- 
tocorrelation function is very small and noisy. To overcome this obstacle, we then fit 
the autocorrelation function to the transient behavior at intermediate times, which will 
overestimate the gap, thereby providing an upper bound. 

As shown in Fig. ID. 51 the gap obtained from the asymptotic behavior for small 
system sizes agrees perfectly with the exact diagonalization results. For intermediate 
system size the transient behavior overestimates the gap. However, for very large sys- 
tem sizes we see that this upper bound goes to zero with increasing system size faster 
than the gap extrapolated from the exact diagonalization results. Fitting the large-A^ be- 
havior to a power-law ' ' we obtain z - 0.765(6) for the loop gas, and z - 0.746(4) 
for the string net. The graph Laplacian times the system size {NX.') is hence gapless. 
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a) loop gas 




Figure D.4: Autocorrelation functions: Monte Carlo results for the autoconelation function of a) the tree 
height for the loop gas b) graph diameter for the string nets. The Monte Carlo time r is given in arbitrary 
units. Lines denote exact diagonalization results. 



There is a simple heuristic argument for the crossover scale between the gapped 
behavior for small and the gapless behavior for large in the case of the string net. 
This argument is best discussed in the dual picture of triangulations of the sphere. For 
a small number of triangles the geometry is always that of a simple sphere and the 
updates mix well resulting in gapped behavior For larger than about 40, one can 
- for the first time - find triangulations that correspond to a geometry of two spheres 
described by two icosahedra connected by a narrow neck. Updates no longer mix well, 
in particular there is a slow mode associated with shifting triangles from one sphere to 
the other via the narrow neck. It is this slow mode which dominates the mixing times 
for large system sizes resulting in gapless behavior. 

Outer planar triangulations 

Finally, in Fig. ID.6l we show Monte Carlo results for the gap of the graph Lapla- 
cian for outer planar triangulations as a function of the number of triangles A^. The 
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Figure D.5: The gap of the graph Laplacian as a function of tlie inverse system size i/N obtained from 
exact diagonalization and Monte Carlo simulations. Results for the off-lattice loop gas are shown on the left 
(1/c' = 3200A'), and for the olf-lattice string net on the right (1/c' = 30A'). The finite-size extrapolation for 
large system sizes reveals that the gap closes as A?"'-^*'**) for the loop gas and as A?"l '''*W for the string net. 
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Figure D.6: The gap of the graph Laplacian for outer planar triangulations (1/c' = N). 



observed decrease as ^ is consistent with but faster than the bound derived in 



Appendix B This scaUng is fast enough to marginally show the gaplessness of any 
local Hamiltonian based on F moves for this model. 
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